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Nitrogen-vacancy (NV) centers in diamond have emerged as valuable tools for sensing and polar¬ 
izing spins. Motivated by potential applications in chemistry, biology and medicine we show that 
NV-based sensors are capable of detecting single spin targets even if they undergo diffusive motion 
in an ambient thermal environment. Focusing on experimentally relevant diffusion regimes we derive 
an effective model for the NV-target interaction, where parameters entering the model are obtained 
from numerical simulations of the target motion. The practicality of our approach is demonstrated 
by analyzing two realistic experimental scenarios: (i) time-resolved sensing of a fluorine nuclear 
spin bound to an N-heterocyclic carbene-ruthenium (NHC-Ru) catalyst that is immobilized on the 
diamond surface and (ii) detection of an electron spin label by an NV center in a nanodiamond, 
both attached to a vibrating chemokine receptor in thermal motion. We find in particular that the 
detachment of a fluorine target from the NHC-Ru carrier molecule can be monitored with a time 
resolution of a few seconds. 

PACS numbers: 76.70.Hb, 75.75.Lf, 76.60.Lz 


I. INTRODUCTION 

Single-molecule observation techniques have provided 
many important insights into the details of fundamental 
chemical and biological processes. Specifically, the in¬ 
vestigation of single-molecule catalytic reactions in var¬ 
ious environments has been made possible by several, 
primarily spectroscopic methods. To give some exam¬ 
ples, fluorescence microscopy has been used to observe 
single-molecule enzymatic dynamics mm and to deter¬ 
mine the spatial distribution of catalytic activity on crys¬ 
tal surfaces [4] . Moreover, scanning tunneling microscopy 
has demonstrated its potential for monitoring oxidation 
catalysis in real time at liquid-solid interfaces [5]. 

More recently, nitrogen-vacancy (NV) centers in bulk 
diamond [5] have been used for spin resonance spec¬ 
troscopy with a sensitivity ultimately reaching the single¬ 
molecule level which makes them a promising new 

tool for probing the dynamics of chemical processes. This 
is particularly true for single-site catalytic reactions m, 
in which case the catalyst is immobilized on a surface 
while maintaining its full catalytic functionality. The 
probing of catalysis on a diamond surface can be accom¬ 
plished, in principle, by using surface-implanted NV cen¬ 
ters to monitor single nuclear spins of the catalyst (see 
Fig. 0 , which yields valuable information about the cat¬ 
alytic reaction, e.g., about conformational changes. As 
the probing method is label-free it is also suitable for 
more general biosensing applications mm- While this 
prospect is certainly appealing, the important question 
remains of whether NV-based sensing of nuclear spins is 
still efficiently achievable under ambient thermal condi¬ 
tions, i.e., when the nuclear target spins exhibit random 
thermal motion. 

In this paper we show that the detection of nuclear 
spins (and target spins in general) undergoing thermal 
diffusion is indeed feasible under realistic experimental 
conditions. We focus our attention on the nuclear spin of 


a fluorine atom covalently bound to a ruthenium-based 
catalyst (NHC-Ru complex [15]) and analyze the time- 
resolved sensing of a possible detachment of the fluo¬ 
rine from the catalyst. Our theoretical study is strongly 
motivated by recent experimental progress: Decoupling 
schemes for the electron spin of the NV now offer the 
possibility to resonantly amplify interactions between the 
NV center and selected spin targets and to suppress deco¬ 
herence caused by surrounding spins. In this way, single 
nuclear spins located in the diamond lattice are routinely 
detected by employing either pulsed or continuous dyna- 



FIG. 1. Sensing of a single fluorine atom undergoing random 
thermal motion using a single near-surface nitrogen-vacancy 
center in bulk diamond: The fluorine-containing molecule 
(N-heterocyclic carbene-ruthenium complex) is covalently at¬ 
tached to the diamond and rotates freely around its an¬ 
chor point, confining the thermal motion of the fluorine to 
a spherical surface. The nuclear spin of the fluorine and the 
microwave-driven electron spin of the NV center interact res¬ 
onantly at the Hartmann-Hahn matching condition. Time- 
resolved sensing of the fluorine in thermal motion is achieved 
by observing the polarization transfer between the spins. 
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mical decoupling HMg. The achieved sensitivity of NV 
centers has moreover been shown to be sufficient to dis¬ 
cern even single spin targets or molecules on the diamond 
surface !ZHg and also to detect spin ensembles surround¬ 
ing NV centers hosted in nanodiamonds m- 

To make the analysis more definite we consider contin¬ 
uous dynamical decoupling of the NV center in combina¬ 
tion with the Hartmann-Hahn double resonance (HHDR) 
scheme [[18], which has been successfully used for detect¬ 
ing single static target spins m- At first we substan¬ 
tially extend the existing theory on the HHDR scheme 
for NV centers moon by including the effects caused by 
the motion of the target spin. Making use of a stochastic 
description of the diffusive target we derive an effective 
quantum master equation for the evolution of the spin 
system formed by the NV center and the target. We find 
in particular that the evolution of the spin system for a 
target undergoing fast diffusion is remarkably similar to 
the case of a static target, but with an effective averaged 
NV-target coupling replacing the static NV-target cou¬ 
pling. Details of the motion of the target are obtained 
from anisotropic network model (ANM) simulations [201 - 
[22] of the entire NHC-Ru complex. The results of the 
ANM simulation allow us to generate stochastic trajec¬ 
tories of the target spin only, which in turn can be used 
for estimating the parameters of the effective quantum 
master equation. 

When applied to the fluorine attached to the NHC-Ru 
complex our approach yields that the effective NV-target 
coupling can reach values up to a few hundred kilohertz 
and therefore is sufficiently strong for observable HHDR 
polarization transfers between the spins to be achieved. 
The motion of the target spin also results in diffusion- 
induced decoherence of the NV spin, which is however 
negligible. As a result, the detachment of a target spin 
from the NHC-Ru carrier molecule can be monitored with 
a time resolution of a few seconds. In addition, we ana¬ 
lyze an alternative scenario, where an electron spin label 
is attached to a CXCR4 chemokine receptor, which plays 
an important role in infection processes |23| . Unlike in 
the previous application, the NV center used to detect 
the spin label is inside a nanodiamond that is attached 
to the same receptor. Our results indicate that spin la¬ 
bels at a distance of 10 nm from the NV center can be 
monitored with a sub-second time resolution despite the 
thermal vibrations of the chemokine receptor. 

The paper is organized as follows: In Sec. |TT] a detailed 
description of the system and the Hartmann-Hahn reso¬ 
nance scheme is presented. In Sec. |III| we characterize the 
motion of the NHC-Ru complex by means of ANM sim¬ 
ulations and introduce stochastic differential equations 
for generating stochastic trajectories of the target spin. 
In Sec. IIVI we start from the stochastic Hamiltonian of 
the NV-target system and derive the quantum master 
equation for the coupled spins. In Sec. [V] we discuss the 
feasibility of the spin detection considering the relaxation 
time of the NV spin and the NV readout photon collec¬ 
tion efficiency. In Sec. |VI| we briefly analyze the detection 
of an electron spin label on a CXCR4 chemokine receptor 


by using a proximate NV center in a nanodiamond. We 
end with the conclusions in Sec. m 

II. MODEL AND STOCHASTIC APPROACH 

Starting with our main application, we consider a fluo¬ 
rine atom bound to the NHC-Ru molecule, which is cova¬ 
lently attached to the surface of the diamond and rotates 
freely around its anchor point (see Fig. [lj. In general, 
the carrier molecule can be flexible, but the NHC-Ru 
complex considered here exhibits only small changes in 
shape. The NHC-Ru molecule is immersed in a layer of 
water, deposited uniformly on the surface of the diamond, 
and undergoes fast diffusive motion due to thermal fluc¬ 
tuations (the environment being at room temperature). 
The NV center is implanted close to the diamond surface 
and in the vicinity of (but not necessarily directly under¬ 
neath) the anchor point. We do not expect to have full 
control over the positions so that the details of the setup 
rely on the statistical proximity between the NV center 
and the molecule. The axis of the NV center is assumed 
to be perpendicular to the diamond surface for concrete 
applications in Sec. [V] 

A. Coupling between NV center and target spin 

The NV center and the fluorine atom (i.e. the target 
spin) are coupled through the magnetic dipole-dipole in¬ 
teraction between their electronic and nuclear spins, re¬ 
spectively. The NV ground state is an electronic spin 
triplet (spin S = 1) with three states having spin projec¬ 
tions | m s = 0) and \m s = ±1) along the NV axis, defin¬ 
ing the z-axis of our coordinate system. In what follows, 
the states \m s =0) and |m s = ±1) will be denoted as |0) 
and |±1), respectively. The target spin is a nuclear spin 
doublet (spin S = 1/2) with spin projections |f) and |j,). 
The degeneracy of the states | — 1), |+1), both separated 
from |0) by the zero-field splitting A = 2.87 GHz, is lifted 
by the external magnetic field B such that a continuous 
microwave field resonant with the |0) -O- | —1) transition 
can be applied (see Fig. [I]). The NV center is then de¬ 
scribed within the 10), | —1) subspace by the microwave 
dressed states |±) = (|0) ± |—l))/v^2 separated by twice 
the Rabi frequency fl [24] . 

The time-varying Hamiltonian used to describe the 
coupled spin system in the secular approximation and 
the dressed state basis is given by |16j [19] 

H(t) = flat + 7 N-B c ff(t) • er N - a%[A(t) • er N ], (1) 

where A(t) is the hyperfine vector, B e s(t) = B — 
| 7 I / 1 A(t) is the effective magnetic field and 7 n is the 
nuclear gyromagnetic ratio. The spin operators cr® and 
a^ (with j = x,y,z) act on the states |+), |—) and |f), 
[/.), respectively. The time dependence of the hyper¬ 
fine vector A(t) stems from the fluctuating position of 
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FIG. 2. Energy-level diagram of the electronic spin 10), | — 1) 
of the NV center and the nuclear spin |d), 4) of the fluorine 
atom, (a) The mismatch between the electronic and nuclear 
energy scales suppresses the polarization transfer between the 
states 10,4.) and | — 1,')') induced by the hyperfine coupling 
A. (b) Continuous microwave driving of the electronic states 
enables polarization transfer between the dressed states |+,4_) 
and | — ,t) under the Hartmann-Hahn matching condition. 


the fluorine r(t) relative to the NV center. Specifically, 
A(t) = A[r(i)] and r(t) are related at every instant of 
time by 


A{r) 
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1 ), 


( 2 ) 


where fj are the components of the unit vector r = r/|r|. 
Moreover, q e and 7 n are the electronic and nuclear gyro- 
magnetic ratios, respectively, and /io is the vacuum per¬ 
meability. 

The Hamiltonian in Eq. (|T]) has been shown to accu¬ 
rately describe the standard HHDR scheme ns, which 
has been employed in a recent experiment to detect 
static target spins by using a single NV center nu 
Specifically, under the Hartman-Hahn matching condi¬ 
tion fl — 7 N |-B cff | = 0, the states |+,i) and | — ,t) of the 
coupled spin system become resonant, which results in 
coherent polarization transfer between them (see Fig. [2]) . 
In the ideal HHDR scheme, the NV spin is first optically 
pumped into the state |0) and rotated to |+) with a 7t/2 
microwave pulse. If the target spin is in state 4) then 
the system is thus prepared in state +, J,) and undergoes 
coherent oscillations between the |+4) and | —,t) states. 
The polarization transfer resulting in the state |—, f) af¬ 
ter the interrogation time r; nt serves as experimental in¬ 
dicator for the presence of the target spin, observed as 
modulations in the NV fluorescence after mapping the 
states |+), |—) back to |0), | —1) with a second n/2 pulse. 
The preparation, polarization transfer and optical read¬ 
out of the NV spin are in principle identical for static and 
diffusive target spins. In practice, however, the transfer 
is affected by different decoherence mechanisms and the 
target spin is initially in a mixed state, which will be 
shown to reduce the observed polarization transfer. 


This motivates our main assumption, namely that the 
position r(t) and the hyperfine vector A(t) are stochas¬ 
tic processes, i.e., random functions of time. As a conse¬ 
quence, the time evolution governed by H(t) becomes a 
stochastic differential equation. 

To make the problem more tractable we restrict our¬ 
selves to the case where r(t) and A(t) are strictly station¬ 
ary stochastic processes [55]. This implies that A(t) has a 
stationary probability distribution pa(s) with mean A = 
(A(t)) p and variance er 2 = (( A(t ) — -4.) 2 ) p , where (-)p 
denotes averaging with respect to pa{ s )- We further as¬ 
sume that A{t) is fully specified by the two-point correla¬ 
tions Cij(r) = (£i(t + T)£j(t))tr expressed in terms of the 
fluctuations of the hyperfine vector (t) = Aj (t) — Aj , 
where now the average (■ ) tr is taken over all diffusive tra¬ 
jectories of the target spin. Instead of using correlations 
Cij(r) we can characterize £(t) in an equivalent way by 
the power spectra 


Sij(u)= / d TCij{r)e lulT . (3) 

J — OO 

For the common case of exponentially decaying correla¬ 
tions Cij(r) = a^je~' T " Ti:i the power spectra are S^( u>) = 
2 a ij T ij /{1 +w 2 Tj 2 ), where afj are the amplitudes of the 
fluctuations and the correlation times. 

Our general approach is to divide the task of deter¬ 
mining the evolution of the coupled spin system into two 
independent parts: The first part deals with the char¬ 
acterization of the diffusive motion of the target spin 
in terms of the distribution pa(s) and the correlations 
Cij(r). The second part consists of solving, at least ap¬ 
proximately, the evolution of the spin system according 
to H(t) with the distribution pa{ s ) and the power spec¬ 
tra Sij(ui) as given ingredients. 


III. DIFFUSION OF THE TARGET SPIN 

To efficiently characterize the diffusive motion of the 
molecule and the fluorine on different time scales we use 
short-time ANM simulations in combination with an ef¬ 
fective description based on stochastic trajectories. More 
precisely, ANM simulations are used to simulate the dy¬ 
namics of the entire carrier molecule, whereas stochastic 
trajectories provide a simple but still accurate descrip¬ 
tion of the random motion of the fluorine or the target 
spin in general. The time T exp required for the fluorine to 
explore the available spatial domain (see Fig. § roughly 
defines the time scale separating the two descriptions. 


B. Stochastic description of diffusive target spins 

When applying the HHDR scheme to the diffusive tar¬ 
get spin we face the additional complication that its posi¬ 
tion r(t) undergoes random fluctuations, i.e., Brownian 
motion, as a consequence of the thermal environment. 


A. Anisotropic network model simulations 

Numerical simulations allow us to obtain physically ac¬ 
curate details of the motion of the carrier molecule to 
which the target spin is bound, either the NHC-Ru com¬ 
plex or the chemokine receptor. Apart from ANM simu¬ 
lations, several alternative numerical methods have been 
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proposed to predict molecular motion, including molec¬ 
ular dynamics (MD) [26] ;57j and normal mode analysis 
(NMA) [25]. ANM simulations are particularly suitable 
for our purposes for the following two reasons: Consider¬ 
ing the large disparity in time scales between the molec¬ 
ular motion and spin evolution we are interested only in 
the slowest molecular degrees of freedom, i.e., rotations 
and low-frequency vibrations, which have been shown to 
be accurately reproduced by the ANM. Moreover, ANM 
simulations are computationally less costly than fully mi¬ 
croscopic MD simulations and therefore applicable to bi¬ 
ologically relevant molecules of considerable size. 

In essence, the ANM is a coarse-grained description of 
the collective dynamics of the molecule, independent of 
the atomic details [2UH22) . Regardless of specific atomic 
interactions, all atoms within a chosen cutoff distance R c 
interact via pairwise harmonic potentials whose strength 
is characterized by a phenomenological spring constant k, 
assumed to be the same for all atom pairs [2D:. The values 
of R c and k are obtained by comparing ANM results to 
MD simulations and experimental data [21| , where larger 
cutoffs R c (implying more interacting pairs) can be com¬ 
pensated by smaller spring constants k [22). The only re¬ 
maining parameters entering the simulation are the equi¬ 
librium positions of the atoms forming the molecule, de¬ 
picted in Fig. [T] for the NHC-Ru complex. The collisions 
of the molecule with the constituents of the solvent at am¬ 
bient temperature are mimicked by augmenting the ANM 
simulations with a Langevin dynamics, i.e., damped ran¬ 
dom kicks, for each atom of the molecule [2D]. In this 
way, we have complete access to the trajectories of the 
target spin, which depend on the shape and stiffness of 
the molecule, the solvent and its temperature [5DU52] . 

When applied to the NHC-Ru molecule we use the val¬ 
ues for the spring constant n = 1 kcal/(mol A 2 ), the cut¬ 
off distance R c = 10A, the damping coefficient of the 
Langevin dynamics </ = 5ps -1 and the solvent temper¬ 
ature T = 300 K. The ANM simulations with this set 
of parameters have been tested against MD simulations 
by comparing vibration spectra of various molecules and 
are expected to give reliable results for the rotations of 
the NHC-Ru complex and the low-frequency vibrations 
of CXCR4 chemokine receptor. The molecular geome¬ 
try and the atomic numbers of the constituents of the 
NHC-Ru complex are provided in XYZ file format as 
Supplemental Material 331 . 

Figure [3] shows representative sample trajectories of 
the fluorine atom attached to the NHC-Ru complex on 
the time scale of microseconds. The fluorine undergoes 
a diffusive motion within a narrow spherical shell whose 
thickness is determined by small oscillations of the NHC- 
Ru molecule. The motion is restricted to a spherical 
zone as a result of the presence of the diamond surface 
and the spatial extent of the molecule. If we neglect the 
thickness of the shell then the position of the fluorine is 
conveniently described in a spherical coordinate system 
centered at the anchor point of the molecule, with az¬ 
imuthal angle ip £ [0, 27t] and polar angle 9 £ [0, n/2]. 
Within this parametrization, we extract the properties 




FIG. 3. Trajectories of the target spin obtained from short- 
time ANM simulations, (a) The fluorine atom bound to the 
NHC-Ru molecule follows a stochastic trajectory on a spheri¬ 
cal surface, (b) The spin label attached to the chemokine re¬ 
ceptor CXCR4 exhibits damped fluctuations around its equi¬ 
librium position with respect to the NY center. 


of the diffusive motion in the form of the distributions 
pe(s) and p v (s) of the polar and azimuthal angle, re¬ 
spectively, and the rotational diffusion coefficient D r , as 
shown in Fig. [4] In particular, the mean square displace¬ 
ment (A9 2 (t))ti in the regime t <C T exp allows us to deter¬ 
mine the diffusion coefficient T> r from the short-time re¬ 
lation (A0 2 (f)) tr = 2V r t. The value of the diffusion coef¬ 
ficient extracted from ANM simulations is V v = 2.1ns _1 
for the specific NHC-Ru molecule immersed in water at 
room temperature. 

We can also estimate the diffusion coefficient 2? r by re¬ 
sorting to semi-empirical formulas, either as a practical 
alternative to ANM simulations or to gain some intu¬ 
ition about the molecular motion. If the approximately 
rod-shaped NHC-Ru molecule is modeled as a cylinder of 
length L and diameter d and the attachment to the di¬ 
amond surface is neglected then the rotational diffusion 
coefficient may be obtained from 130 

Pr = ^^ [l0gP + C(:p)] ’ (4) 

with c(p) = —0.05/p 2 + 0.917/p — 0.662. The diffusion 
coefficient T> r depends on the ratio p = L/d , the Boltz¬ 
mann factor h-gT and the viscosity of the solvent rj. The 
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coefficient T> r depends crucially on the length L , imply¬ 
ing that larger molecules undergo slower diffusion. For 
a molecule of size L = 1.5 nrn and d = 1.1 nm dissolved 
in water at T = 300 K with viscosity 77 = 10~ 3 Pas we 
obtain T> r = 0.3 ns -1 , thereby underestimating the more 
accurate simulation result. 


B. Stochastic trajectories and correlation functions 


In order to specify the correlations Cjj(r) we require 
stochastic trajectories on longer time scales than acces¬ 
sible by numerically costly ANM simulations. For that 
purpose we generate stochastic trajectories of the fluo¬ 
rine on a spherical surface using the stochastic differential 
equations in Ito form |34| : 


d 0 (f) 
d <P(t) 


Vr 

tan 0 ( 1 ) 


dt + v^dWitf) 


y/2V~ r 
sin 0 (f) 


d W 2 (t), 


(5) 


where dWi(t) and dW 2 (t) are two independent Wiener 
processes and dt is the time differential. In this setting, 
we neglect the spatial extent of the molecule and restrict 
the motion of the fluorine to a spherical zone limited by 
the polar angles 0 m ; n and 0 max , which is implemented 
through reflective boundary conditions for Eqs. ©■. The 
distribution of the fluorine on the spherical zone resulting 
from Eqs. (|5j) is uniform; the corresponding distributions 
of the angles are p v {s) = l/2n, with s € [0, 27t], and 

. sin(s) 

P9 S cos( 0 mi „) - cos( 0 max ) ’ 

with s G [0min,0max]- The angles 0 min and 0 max are cho¬ 
sen to best fit the ANM results for pg(s), as illustrated in 
Fig. gb). However, the results of our analysis are rather 
robust with respect to the choice of 0 m i n and 9 max , which 
alternatively may be estimated based on the shape of the 
molecule. 

Now, for a generic stochastic process u(t) evolving ac¬ 
cording to d u(t) = a[it(f)]df + b[u(t)]dW(t), with drift 
coefficient a[u(t)] and diffusion coefficient b[u(t)}, one ob¬ 
tains numerically generated trajectories by using the cor¬ 
responding finite difference equation 

u{t + St) = u{t) + a[u(t)]5t + 6[u(t)](5t) 1 ^ 2 7V(t), (7) 


where St is the discrete time step and N(t) is a tempo¬ 
rally uncorrelated normal random variable, i.e., N(t) is 
statistically independent of N(t'). We apply the finite 
difference approach to generate trajectories [ 0 (f), tp(f)] 
based on Eqs. © with the parameters 0 m j n , 0 max and 
T> r obtained from ANM simulations. Trajectories of the 
position r(t) relative to the NV center and the hyperfine 
vector A(t) are then readily found from [0(f), <p(t)\. 

The two-point correlations CV,-(t) (with i,j = x,y,z ) 
obtained from Monte Carlo sampling over trajectories are 


shown in Fig. [5j While the details of C'jj(r) depend on 
the random positioning of the NV center, we find that 
the correlations decay quickly with time. This ensures 
that the power spectra Sij(ui ) and integrals over the two- 
point correlations Cij(r) are well defined, and that the 
Markov approximation is applicable in Sec. |IVB| In fact, 
the specific dependence of the stochastic equations © 
on the diffusion coefficient V r implies that all correla¬ 
tion functions CV,(t ) decay on a time scale set by V~ x 
and thus ~ V~ Y for exponentially decaying (r) 
(see the Appendix for details). 


IV. EVOLUTION OF THE COUPLED SPINS 


For the analysis of the stochastic evolution we split 
the Hamiltonian H(t ) into a time-independent average 
Hamiltonian H= (H(t)) p and the remaining stochastic 
part H^(t) = H(t) — ( H(t )) p , which read as 


Ha = nB - • er N - <j e x [A ■ a 


w) = 


-o-xK(i) 


( 8 ) 


We note that Ha is identical to the Hamiltonian for a 
static target spin, where the constant hyperfine vector 
A is replaced by the mean A = (A(t)) p , while H^(t) 
depends only on the fluctuations £(t) of the hyperfine 
vector. The evolution of the density matrix p{t) of the 
spin system is governed by the Liouville equation 


^_p(t) = [C A + Cz(t)]p(t) 


(9) 


where the Liouvillian operators are defined as C A = 
— i [Ha, • } and C$(t) = —i [H^(t), ■ ]. As it stands, Eq. ([ 9 ]) 
describes the evolution of the coupled spins for a sin¬ 
gle diffusive trajectory of the target spin, encoded in the 
fluctuations £(i). Since we can control neither the ini¬ 
tial position nor the trajectories of the target spin our 
goal is to determine (p(i))tr, i- e -j th e density matrix p{t) 
averaged over all possible trajectories £(i) with random 
initial conditions. 

A straightforward way to calculate (p(t)) t r approxi¬ 
mately is by Monte Carlo sampling: For a representative 
ensemble of trajectories {£(i)| one solves Eq. © numer¬ 
ically for each £(i) and takes the ensemble average of the 
resulting density matrices to obtain (p(f)) tr . The trajec¬ 
tories £(t) may be generated by following the procedure 
outlined in Sec. Ill B We will use Monte Carlo sampling 
to confirm and complement our analytical results at the 
end of this section. 

Approximate analytical descriptions for the evolution 
of {p(t ))tr can be obtained for several relevant parameter 
regimes of the system. The starting point is the stochas¬ 
tic Liouville equation [551 125] 




f dr(C^(t)e CAT C^(t — r)) tr e 
Jo 


-Cat 


(10) 


(p(t))t 
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FIG. 4. The properties of the diffusive motion of the fluorine (target spin) extracted from the ANM simulations, (a) The 
mean square displacement ( A6 2 (t))tr as a function of time t yields the diffusion coefficient T> r = 2.1ns _1 from the short-time 
relation {A9 2 (t))t T = 2V T t. The histograms (b) and (c) show the distributions pg(s ) and Pip(s) of the polar and azimuthal angle, 
respectively. The dotted lines in (b) and (c) indicate p v (s) and pg (s) expected from a uniform distribution of the fluorine on a 
spherical zone limited by the polar angles 9 m i n and # max . 


which is obtained by formally integrating the Liouville 
equation ([9]) and performing a cumulant expansion up 
to second order in C^{t). The applicability of Eq. (10) 
is limited by the amplitude of the fluctuations £(f) and 
their correlation times, which can be quantified by the 
typical amplitude and width of CV,(t), denoted by <f 2 
and f, respectively. The cumulant expansion includes 


terms up to order <r 2 f f, which implies that Eq. (10) is 


valid under the condition <f 2 f <C T; -1 for spin evolutions 
up to the interrogation time T; nt . 

We now discuss three different regimes for which an¬ 
alytical approximations are obtained: the fast-diffusion 
regime with f —> 0 and arbitrarily large a 2 , the moderate- 
diffusion regime with finite a 2 and f, and the slow- 
diffusion regime with f Tint- 


where the Pauli matrices &j act on the subspace only. 
For the nuclear spin of the fluorine in a completely mixed 
state and the NV center initially in the state |+), we then 
obtain the probability P(f, 6) of finding the NV center in 
the reversed state |—) after the time t 

Pit, 6) = \j^T¥ s ™ 2 ( 12 ) 

where we introduced the averaged coupling between the 
NV and the target spin J = j{A 2 + A 2 ) 1 / 2 . The maxi¬ 
mal polarization transfer occurs at resonance 5 = 0 after 
a time t = 7r/2JP, which implies that the interrogation 
time T; nt is of the order of 1/ J , typically in the millisec¬ 
ond regime (see also Sec. [V|. 


A. Fast-diffusion regime 


B. Moderate-diffusion regime 


For vanishingly small correlation times t —> 0 the sec¬ 
ond order term in Eq. (10) and higher order terms of the 
cumulant expansion vanish. The evolution of the spin 
system is governed solely by the average Hamiltonian Hy[ 
in Eq. (J8| . The only difference from a static target spin is 
that the hyperfine vector A is replaced by the mean A of 
the stationary process A{t) over the distribution pa(s). 
This simplifies the description of the system considerably 
because no averaging over trajectories is required. 

Analogously to the static case mm we determine 
the dynamics of the spins according to Hy i (with the 
magnetic field B along the NV axis). In the vicinity 
of the Hartmann-Hahn matching condition, i.e., close to 
vanishing averaged detuning S = ^(f2 — 7 n-B z + ^A z ), 
transitions to the off-resonant states |+,t) and | — ,i) are 
suppressed by the factor (A/fl) 2 of the order 10 - °. This 
allows us to reduce the system to the two-dimensional 
subspace |+,4_), | —, t) with the corresponding Hamilto¬ 
nian 


V.4 — 2 do z ^ A x Ax 


2 AyOy , 


In the moderate-diffusion regime with a finite correla¬ 
tion time f, one has to take the second order correction 
in Eq. (10) into account. However, since the evolution 


according to Hy± is slow on a time scale set by f, i.e., 
Jt <C 1, we can neglect the exponentials e ±CAT and ob¬ 
tain 


6' 
dt 


{pit)) tr = (ca + J o d T{£ ( {t)C ( {t - T)) tr ^ 


{pit)) t 


(13) 

The effect of the stochastic part is best understood when 


Eq. (131 is reduced to a Markovian master equation in 
Lindblad form [Ml EH- In the Markov approximation 
we extend the upper limit of the integral in Eq. (13) to 


infinity and introduce the rates yij = Sij(uj = 0), which 
for the particular case of exponentially decaying corre¬ 
lations Cij(r) = <7 2 -e - l T l/ r o are given by 7 ^ = 2er 


Tin • 


ij 1 IJ 


More generally, we note that the matrix T constituted by 
the elements 7 ^ is positive definite for any form of the 
correlations C i: j (t) [55]. After a change to the eigenbasis 
( 11 ) of the symmetric matrix T, having eigenvalues 7 k, the 
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( a ) 



r(ns) 


|—), which directly affect the Hartmann-Hahn polariza¬ 
tion transfer. In fact, the dissipative part results from 
averaging over dissimilar trajectories £(t), which during 
the interrogation time T\ nt explore slightly different do¬ 
mains accessible to the hyperfine vector. This is a conse¬ 
quence of the finite correlations in combination with the 
random initial values of £(t). On the other hand, the po¬ 
larization transfer when only a single diffusive trajectory 
of the fluorine is considered remains always coherent. 


C. Slow-diffusion regime 


(b) 



r(ns) 

FIG. 5. Characterization of the thermal fluctuations £(t) of 
the hyperfine vector. The correlations Cij(r) are normalized 
by the standard deviations <t; of the fluctuations, (a) The au¬ 
tocorrelations Cu(t) are, to a good approximation, exponen¬ 
tial decay curves with a typical correlation time of t,^ ~ 0.1 ns. 
(b) The cross-correlations Cij{r) with i A j are found to be 
significantly smaller than the autocorrelations. 


resulting Lindblad equation reads as 

d 


dt 


(p(t))tr = -i [Ha, {pit)) tr] + £diff(p(X))tr (14) 


with the diffusion-induced dissipative part 

-£diff {pit)) Ir = (^Lk{p{t))tiL\~ {p{t)) tr}^ , 

k 

(i5) 

where { , } stands here for the anticommutator. The three 
Lindblad operators in Eq. (|15|) acting on both spins are 


given by Lk = L' k = <g» (er® + Ifb), where dj) 1 are the 

Pauli matrices along the axes defined by the eigenbasis of 
T. This effective description of the coupled spin system 
in the form of Eqs. (14) and (151 greatly simplifies the 


analysis of the HHDR scheme and is applicable for a large 
range of the detuning S. 


It can be seen from Eq. (14) that the coherent evolu¬ 


tion is determined by the averaged Hamiltonian H Al as 
in the fast-diffusion regime, while the diffusion-induced 
dissipative part of Eq. (14) has two effects: The opera¬ 
tors act in three spatial directions on the nuclear spin 
and result in its depolarization. The operator cr® act¬ 
ing on the NV center causes spin flips between |+) and 


In the slow-diffusion regime we assume that the hy¬ 
perfine vector A{t) is constant during the interrogation 
time ri n t, but drifts slowly during repeated measurements 
of the Hartmann-Hahn polarization transfer. This means 
that Ti n t -C f -C IWint, where N is the number of individ¬ 
ual measurement runs that are performed. This regime 
was also analyzed in the context of magnetometry with 
rotationally diffusing nanodiamonds j59| . 

To obtain the averaged density matrix {p(t )) tr we 
determine the spin dynamics according to the time- 
independent Hamiltonian H as for the static case, but 
with the hyperfine vector A as a random variable with 
distribution pa( s )- This yields the density matrix p(t, A) 
for each A , from which we obtain 

(p(*))tr = j dsp A (s)p(t,s) . ( 16 ) 


The density matrix (p(t)) tr is expected to accurately de¬ 
scribe the observed experimental results in the limit of 
large N. 

To emphasize the difference from the fast-diffusion 
regime we give the explicit expression for the probability 
Pit, 6) of the transition from state |+) to |—). According 
to the prescription in Eq. (16) this probability is given 

by 


p {t,ty = l I dsp A (s) J sin 2 [v{s)t], 

* J F(s) + d 2 (s) 

( 17 ) 

where the detuning 5{A) = 1(H — 7 n-B z + \A z ), the cou¬ 
pling J{A) = + Ay j 1 / 2 and the effective Rabi fre¬ 

quency v{A) = [J 2 {A) + 5 2 {A)] 1 t 2 are random variables. 
In the slow-diffusion regime we thus solve the equation of 
motion and subsequently average the solution over p A is), 
in contrast to the fast- and moderate-diffusion regimes, 
where the equation of motion is solved for averaged de¬ 
tuning S and coupling J. 

Figure [G] shows the polarization transfer probabili¬ 
ties P(Tj nt ,fl) for the different diffusion regimes. The 
general agreement between the analytical results and 
Monte Carlo sampling over trajectories is illustrated in 
Figs. |6](a) and[6](b) for a set of generic parameters. The 
analytical results in Fig. [6](a) are accurate in the regime 
Tjnt lms, as expected from the condition <t 2 t <C tT* 




















H(MHz) 

FIG. 6. Polarization transfer probabilities P(ri n t,fi) in the 
fast-, moderate- and slow-diffusion regimes, obtained from 
Monte Carlo simulations (symbols) and analytical results 
(lines), (a) Transfer probability P(n n t,fl) for fixed Rabi fre¬ 
quency SI corresponding to a detuning 5 = 1.8 kHz. Analyti¬ 
cal results and Monte Carlo sampling over trajectories are ex¬ 
pected to agree in the regime Tj n t -C 1 ms. (b) Transfer prob¬ 
ability P(rint, fi) for a fixed interrogation time Tint = 0.25 ms. 
(c) Maximum achievable transfer probability PMax(SI) for 
given Rabi frequency S2 and optimized interrogation time 
Tint- Perfect polarization transfer PMax(SI) = 0.5 is achieved 
for fast diffusion, while moderate and slow diffusion result 
in reduced transfer. V T = (0,10,10 3 ) ps _1 in (a), (b) for 
the slow-, moderate- and fast-diffusion regimes, and V v = 
(0,1,10 3 ) ps' 1 in (c), with random initial positions of the 
target spin. 


for the validity of the second order cumulant expan¬ 
sion. The comparison in Fig. | 6 ])b) demonstrates that the 
analytical approximations are valid even far from reso¬ 
nance. Finally, the maximum achievable transfer proba¬ 
bility -PMax(^) f° r a given Q and optimized interrogation 
times Ti n t is shown in Fig. | 6 |c). Maximal polarization 


transfer is achieved in the fast-diffusion regime, whereas 
reduced polarization occurs in the moderate- and slow- 
diffusion regimes due to diffusion-induced decoherence of 
the NV spin. The most detrimental effect is observed for 
slow diffusion, where the time scale of the diffusion I ?” 1 
and the interrogation time Ti nt are comparable. Note that 
the resonance of Pm & x(^) for fast diffusion (and to a good 
approximation for moderate diffusion) corresponds to a 
Lorentzian of width J centered at <5 = 0, whereas the 
resonance in the slow-diffusion regime is asymmetric and 
not related to the averaged detuning 5. 


V. DETECTION OF FLUORINE 
DETACHMENT FROM THE MOLECULE 


The general formalism developed in the previous sec¬ 
tion is now applied to the concrete problem of detect¬ 
ing the detachment of the fluorine atom from the car¬ 
rier molecule, which demonstrates the time-resolved op¬ 
eration of the HHDR scheme in a thermal environment. 
The situation that we envisage here is that the fluorine- 
marked NHC-Ru complex serves as a catalyst for a chem¬ 
ical reaction. The departure of the fluorine then indicates 
the event of a chemical reaction that led to the exchange 
of the fluorine. 

As previously shown, the diffusive motion of the flu¬ 
orine attached to the molecule (in a solvent at room 
temperature) exhibits short correlation times and the 
Hartmann-Hahn polarization transfer therefore falls into 
the fast-diffusion regime, which is accurately described 
by the corresponding Hamiltonian in Eq. (11). Never¬ 


theless, in order to present the most general situation we 
take the small diffusion-induced decoherence into account 
and accordingly extend the analysis to the moderate- 
diffusion regime. In addition, we include incoherent pro¬ 
cesses caused by spin flips between the states |+) and 
|—), which are the main source of decoherence for shal¬ 
lowly implanted NV centers (40]. The effect of spin flips 
is described by the additional Lindblad dissipator 

£flip<p(t)>tr = ^ E ( L M*)) (p(f)) tr } 

3 


with the Lindblad operators L\ = er®_ and L 2 = cr® , 
which induce spin flips between the states |—) and |+). 
The flip rate 7 m P = 1 /T lp is determined by the relaxation 
time T\p in the rotating frame |41| , being in the millisec¬ 
ond regime for realistic experimental conditions mm- 
The evolution of the spin system is thus governed by the 
Lindblad equation 

^(p(*))tr = -'AHA, (p(*))tr]+>Cdiff(p(t))tr+>Cflip(p(t))tr • 

(18) 

We note that the detrimental effect of spin flips on the 
polarization transfer is significantly stronger than the 
diffusion-induced decoherence given that 7 fl; p 7 *. We 
also recall that the initial state of the fluorine spin is com¬ 
pletely mixed, i.e., p(t = 0 ) = |+)(+|® 5 IF, which reduces 
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(a) 





FIG. 7. The time-resolved detection (with unit signal-to-noise 
ratio) of the detachment of the fluorine from the NHC-Ru 
carrier molecule, (a) The optimal interrogation time Tint for 
each measurement run depending on the flip rate 7 fli p for 
different couplings J. With increasing 7 flj p more runs with 
shorter interrogation times Tint are required, (b) The overall 
measurement time T required to detect the detachment as a 
function of the flip rate 7 fli p for different couplings J. (c) The 
dependence of averaged coupling J on the depth znv of the 
NV center for different lateral offsets xnv with respect to the 
anchor point of the NHC-Ru molecule (see Fig. [I]). 


the observable polarization transfer by approximately a 
factor of two. 

The time-resolved detection based on the HHDR 
scheme is implemented as follows: As long as the flu¬ 
orine is bound to the molecule the probability of finding 
the NV center in state |—) after the interrogation time 
T int is P(Tint,<5). After the detachment the fluorine dif¬ 


fuses away very rapidly from the molecule due to its small 
mass so that effectively J = 0 with the corresponding 
probability Py = o(Ti nt , S). The measurable signal is thus 
A P = P (Ti n t, S) — Pj=o (Rnt, S). The presence of the 
fluorine is confirmed once the signal-to-noise ratio (Psn) 
exceeds unity, which for the optical readout of the NV is 
equivalent to 02] 


A P > 


1 

cVw’ 


(19) 


where N is the total number of measurements and C 
specifies the signal contrast and photon collection effi¬ 
ciency. 

The optimal time resolution is reached by tuning the 
interrogation time T; n t in dependence on the system pa¬ 
rameters, in particular the flip rate 7 fli p . Since the num¬ 
ber of measurements is related to the overall measure¬ 
ment time T through N = T/(r- m t + To), with To <C T- lnt 
the preparation and readout time, we obtain 


T = 


Tint + To 

C 2 A P 2 ’ 


( 20 ) 


from the condition 7 ?sn = 1- The main limitation to the 
time resolution clearly stems from the factor C and the 
flip rate 7 fl ip (see Table [i]) . 

Figure [7] shows the measurement time T and the in¬ 
terrogation time Tint as a function of the flip rate 7 fli p , 


as determined from Eq. (18). The detachment of the 
fluorine can be detected (with unit signal-to-noise ratio) 
with a time resolution of the order of 10 s under realistic 
experimental conditions. The general strategy for mini¬ 
mizing the overall measurement time T in the presence of 
spin flips at rate 7 fli P is to shorten the interrogation time 
T int while increasing the number N of measurement runs. 
More precisely, the interrogation time Tj nt scales as 1 / 7 fli p 
as would be expected from a more generic analysis |44] . 

It can be seen from Fig. [7](b) that short measurement 
times T require sufficiently strong couplings J between 
the NV and the fluorine spin, which in turn are deter¬ 
mined by the depth znv of the NV center and its lateral 


znv 

ZNV 

J 

C* 

7flip + 

Tint 

T 

2.0 nm 

1.5 nm 

0.7 kHz 

0.05 

0.5 kHz 

0.23ms 

2.6s 

2.0 nm 

1.5 nm 

0.7 kHz 

0.05 

1.0 kHz 

0.18 ms 

7.9s 

2.0 nm 

1.5 nm 

0.7 kHz 

0.05 

2.0 kHz 

0.12ms 

37s 

3.0 nm 

1.5 nm 

0.25 kHz 

0.05 

0.5 kHz 

0 - 42 ms 

45s 

3.0 nm 

1.5 nm 

0.25 kHz 

0.05 

1.0 kHz 

0.26ms 

250 s 


* from Refs. PHIS] 
t from Ref. US 


TABLE I. Specific values of the model parameters and mea¬ 
surement times T for detecting the fluorine detachment from 
the NHC-Ru molecule. The coordinates znv and xnv deter¬ 
mine respectively the depth and lateral offset of the NV with 
respect to the anchor point of the molecule (see Fig. [T]) . The 
gyromagnetic ratio of fluorine is 7 f = 40.1 MHz T -1 and the 
preparation and readout time to ~ 100 ns is neglected. 


























10 


offset snv with respect to the NHC-Ru molecule (see 
Fig. 0 As shown in Fig. [7](c) , the coupling J decreases 
monotonically with the depth z^v, where in particular 
an experimentally accessible depth of znv = 2 nm corre¬ 
sponds to J = 0.7 kHz for xnv = 0.5 nm. The depen¬ 
dence of J on the lateral offset Inv with respect to the 
anchor point of the molecule is characterized by the di¬ 
rectionality of the dipolar spin interaction. Qualitatively 
speaking, the coupling J is strong if the averaged posi¬ 
tion of the target spin is located within the dipolar lobes 
of the NV spin, whereas J vanishes for the perfectly sym¬ 
metric configuration ccnv = 0 and for large lateral offsets. 
We finally note that for a fully coherent evolution of the 
spins (7 k = 7 fli P = 0) the optimal interrogation time is 
r int = 7r/2^7, equivalent to a full polarization transfer. 


VI. PROBING VIBRATING MOLECULES 
WITH NANODIAMONDS 


As a second application of our general approach we 
consider an NV center inside a nanodiamond, which is 
attached to a larger molecule, here the chemokine recep¬ 
tor CXCR4 [23]. The aim is to detect an electron spin 
label [S] attached to the same molecule in the vicinity of 
the nanodiamond. There are two main differences from 
the previous scenario: First, the gyromagnetic ratio of 
the electron spin is three orders of magnitude larger than 
that of nuclear spins, which enhances the NV-target cou¬ 
pling considerably. Electron spin labels up to a distance 
of 10 nm from the NV can therefore be detected. Second, 
the relaxation time of the NV spin is shorter in nano¬ 
diamonds than in bulk diamond because of interactions 
with surface spins [39] . 

As for the previous application, we first analyze the 
stochastic properties of the random thermal motion of 
the target spin based on ANM simulations. The parame¬ 


ters entering the simulation are the same as in Sec. Ill A 


and the structural information of chemokine receptor is 
obtained from the Protein Data Bank (PDB) (32- We 
focus on the impact of molecular motion on the relative 
distance r(t) between the NV center and the spin label 
under the assumption that the nanodiamond is station¬ 
ary with respect to the external magnetic field. Figure [3] 
shows representative sample trajectories, obtained from 
ANM simulations, of the relative distance r(t) between 
the NV center and the spin label. The distance r(t) can 
be expressed as r(t ) = R + q(t ), where R is the equi¬ 
librium position and q(t) represents the thermal fluctua¬ 
tions. The numerical results suggest that q(t) can be ac¬ 
curately described by an Ornstein-Uhlenbeck (OU) pro¬ 
cess with the corresponding stochastic differential equa¬ 
tion [25] 


d q k (t) = -rj k q k (t)dt + V2VdW(t), (21) 

where the diffusion is assumed to be isotropic and the 
components q k {t) are defined with respect to a local or¬ 
thogonal basis. The restoring force of strength r\ k ex- 


(a) 




M kHz) 


FIG. 8 . Characterization of the thermal fluctuations £(f) of 
the hyperfine vector, assumed to be an OU process, (a) The 
correlations Cij(r) are exponentially decaying with typical 
correlation times 77 , ~ 1 ps. The movement of the spin la¬ 
bel is approximately isotropic with comparable correlations 
in all spatial directions, (b) The distribution pa x (s) of the x- 
component of the hyperfine vector A(t). The approximately 
Gaussian shape of (s) is a consequence of the underlying 
OU process of the fluctuations q(t). 


erted by the molecule acts against the thermal fluctua¬ 
tions which are characterized by the diffusion coefficient 
V. The stationary distribution of the fluctuations q k [t) 
is Gaussian with variance (A q k (t))ti = R/qk, whereas 
(A^(t)) tr = 2Vt in the regime t <C r exp . We use the two 
relations to determine the parameters q k and V based on 
numerical estimates of (Ag 2 (f)) tr from the simulations. 
Typical values obtained in this way are q k ~ 10 ps -1 and 
T> ~ 10 -5 cm 2 s _1 for an equilibrium distance |_R| = 8nm 
and the environment at room temperature. 

The distribution pa (s) and the correlations CV,(t) are 
obtained from the stochastic trajectories £(i) of the hy¬ 
perfine vector, which are generated by numerically solv¬ 
ing Eq. © with the parameters from ANM simulations 
and subsequently using the explicit expression for A(r) 
in Eq. Figure [8] shows the stationary distribution 
PA x (s) and the correlations CV,(t). The correlations de¬ 
cay on a time scale f ~ 1 ps such that the condition 
<f 2 f <C Tjjjj 1 for the validity of second order cumulant ex¬ 
pansion is easily fulfilled. We note that the correlation 
time is an order of magnitude smaller than in the previ- 















11 


ous example. This is to be expected as the correlations 
Cij(r) resulting from an underlying OU process decay on 
a time scale set by the restoring force rjk rather than by 
the diffusion coefficient V [231 . 

We finally proceed as in Sec. [V] to determine the in¬ 
terrogation time Ti n t and total measurement time T re¬ 
quired to detect the spin label with unit signal-to-noise 
ratio. For the realistic values of the coupling J = 20 kHz, 
the spin flip rate 7 fli p = 1 kHz and C = 0.03 we obtain 
Tint = 0.11 ms and T = 0.6 s. Consequently, the de¬ 
tection of electron spin labels attached to the vibrating 
chemokine receptor in a thermal environment should be 
possible in experiments; with the proviso, however, that 
the rotation of the nanodiamond with respect to the mag¬ 
netic field is negligible [55] , 


VII. CONCLUSIONS 

In this paper we have shown that detection of single 
nuclear and electronic spins undergoing diffusion is feasi¬ 
ble under realistic experimental conditions by using the 
Hartmann-Hahn double resonance scheme. Specifically, 
in the case of rapidly diffusing target spins, which is par¬ 
ticularly relevant for experiments at room temperature, 
we found that the main effect of the diffusion can be ab¬ 
sorbed in an effective NV-target coupling. We moreover 
quantified the decoherence of the NV center induced by 
the target diffusion, which in the limit of fast diffusion 
turns out to be significantly smaller than incoherent ef¬ 
fects caused by surface magnetic noise. 

We emphasize that our approach, namely combining 
a statistical description of the target diffusion with the 
quantum description of the spin dynamics, is very ver¬ 
satile. Instead of using ANM simulations for character¬ 
izing the statistical properties of the target motion one 
may use more sophisticated molecular dynamics simula¬ 
tions [26[ 127) or, alternatively, rely on experimental data. 
This is likely to be necessary for setups involving target 
spins in complex biological environments. 

While the main focus of this paper was on spin sens¬ 
ing, our effective model [cf. Eqs. |l4| ) and ( fl5| )] is also 
perfectly suited for analyzing protocols for the polariza¬ 
tion of external nuclei by means of optically pumped NV 
centers [TH iTSl - EU] . This is particularly true for model¬ 
ing the dynamical polarization of gases and liquids m, 
where the diffusive motion of the target spins plays an 
important role. 
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APPENDIX: RELATION BETWEEN DIFFUSION 
COEFFICIENT AND CORRELATION TIME 


We want to clarify the relation between the correlation 
time of the hyperfme vector A(t) and the rotational dif¬ 
fusion coefficient T> r . For simplicity of the argument we 
consider the scalar case only; the generalization to the 
vectorial case is based on the same reasoning. 

We assume that the hyperfme vector A(t) depends on 
the quantity y(t), corresponding to 9(t) or ip(t) in the 
main text, whose evolution is described by a stochastic 
differential equation of the form 


d y = n(y)dt + a(y)dW. (22) 


The diffusion term a(y) depends on V r as cr(y) ~ \/T\ 
and we assume that the drift term p(y) either depends 
on D r as p(y) ~ T> r or is zero. These assumptions are 
compatible with evolution described by Eqs. ([5]). 

The stochastic differential equation for the hyperfme 
vector A(t) = A[y(t)\ is according to Ito’s lemma 


dA(f) = 


/ dA 1 2 d 2 A\ 

V dy + 2 a By 2 ) 


BA 

dt + a—dW 
dy 


= fi(A)dt + d(A)dW. 


(23) 


Equation (23) has the same form as Eq. (22) with the 


diffusion term d{A) ~ s/TX- The drift term in Eq. (23) 
necessarily depends on T> r as p(A) ~ V r except in the 
trivial case, where y(y) is zero and A(t) is simply pro¬ 
portional to y{t). 

To determine the correlations of A(t) we need the prob¬ 
ability density p(A, t) governed by the Fokker-Planck 
(FP) equation corresponding to Eq. (23), which is 


|p(A,t) = -±[ji(A)p(A,t)] + 1 --^[a\A)p(A,t)]. 

( 24 ) 

The important observation is that the right-hand side of 
Eq. (241 is proportional to T> r such that we can introduce 


the rescaled time t = T> r t. Using the standard methods 
for solving the FP equation (24) with reflecting bound¬ 


ary conditions one then obtains for the autocorrelation 
function EH 


(A(t)A( 0)) tr = £ 


dAAP\(A) 


— Xt 


(25) 


where P\(A) are eigenfunctions with eigenvalues A of 
the time-independent FP equation. Thus, we see from 
Eq. (25) that the correlations of A(t) decay exponentially 


on a time scale set by V r . We note that this result does 
not apply to the OU process in Eq. (21). 


SUPPLEMENTAL MATERIAL 

The XYZ file contains the atomic numbers and the 
Cartesian coordinates of the equilibrium positions of the 
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constituents of the NHC-Ru complex. The atom used to 

1 

4.820 

0.409 

- 2.389 

attach the NHC-Ru complex to the diamond is identified 

1 

3.040 

0.226 

- 2.676 

as a helium atom, corresponding to a carbon atom of the 

1 

3.915 

1.629 

- 3.353 

diamond surface. 



1 

4.651 

5.643 

- 0.001 





1 

3.118 

5.785 

- 0.950 

76 




1 

3.068 

5.834 

0.841 





1 

0.942 

- 1.961 

- 2.679 

6 

0.479 

- 3.582 

1.149 

1 

2.150 

- 3.309 

- 2.682 

6 

- 0.564 

- 4.513 

1.115 

1 

0.561 

- 3.516 

- 3.490 

6 

- 1.087 

- 4.982 

- 0.094 

1 

0.868 

- 1.978 

2.553 

6 

- 0.542 

- 4.505 

- 1.289 

1 

0.510 

- 3.550 

3.341 

6 

0.502 

- 3.575 

- 1.299 

1 

2.104 

- 3.300 

2.557 

6 

1.002 

- 3.130 

- 0.069 

7 

- 1.670 

- 7.323 

- 0.298 

7 

1.994 

- 2.229 

- 0.059 

6 

- 1.912 

- 8.060 

- 1.477 

6 

- 2.214 

- 5.987 

- 0.116 

8 

- 1.410 

- 9.173 

- 1.631 

6 

3.317 

- 2.444 

- 0.046 

2 

- 2.914 

- 7.379 

- 2.740 

6 

3.947 

- 1.202 

- 0.022 

6 

- 1.749 

- 0.740 

0.105 

6 

2.861 

1.787 

1.146 

6 

- 2.592 

0.461 

0.161 

6 

3.011 

3.178 

1.131 

6 

- 1.944 

1.732 

0.106 

6 

3.442 

3.856 

- 0.013 

6 

- 2.736 

2.880 

0.306 

6 

3.719 

3.114 

- 1.165 

6 

- 4.128 

2.782 

0.467 

6 

3.584 

1.722 

- 1.192 

6 

- 4.755 

1.526 

0.463 

6 

3.154 

1.075 

- 0.025 

6 

- 3.984 

0.362 

0.324 

6 

2.387 

1.087 

2.398 

9 

- 4.920 

3.985 

0.651 

6 

3.854 

0.957 

- 2.467 

8 

- 0.516 

1.888 

- 0.041 

7 

2.997 

- 0.255 

- 0.035 

6 

- 0.091 

2.288 

- 1.366 

6 

3.575 

5.358 

- 0.028 

6 

- 0.646 

1.384 

- 2.499 

6 

1.824 

- 0.901 

- 0.051 

6 

- 0.333 

3.790 

- 1.662 

1 

3.749 

- 3.458 

- 0.045 

1 

- 1.073 

- 7.706 

0.404 

1 

5.022 

- 0.958 

- 0.002 

1 

- 2.031 

- 1.743 

0.164 

44 

0.000 

0.000 

- 0.000 

1 

- 2.284 

3.828 

0.341 

6 

1.066 

- 3.066 

- 2.605 

1 

- 5.796 
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1 
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1 
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- 4.886 
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1 
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- 3.447 

1 
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- 2.245 

1 

- 1.726 
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- 2.594 

1 

- 2.796 

- 5.956 

0.834 

1 

- 0.419 

0.340 

- 2.280 

1 

- 2.931 

- 5.732 

- 0.931 

1 

0.006 

4.391 

- 0.817 

1 

2.776 

3.756 

2.039 

1 

- 1.390 

3.982 

- 1.844 

1 

4.041 

3.646 

- 2.074 

1 

0.233 

4.084 

- 2.549 

1 

1.329 

0.761 

2.274 

1 

0.993 

2.150 

- 1.373 

1 

3.013 

0.194 

2.622 





1 

2.434 

1.748 
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